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The correlation function of the two dimensional Ising model with the nearest neigh- 
bours interaction on the finite size lattice with the periodical boundary conditions is 
derived. The expressions similar to the form factor expansion are obtained both for 
the paramagnetic and ferromagnetic regions of coupling parameter. The peculiarities 
caused by finite size are analyzed. The scaling limit of the lattice form factor expansion 
is evaluated. 
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1 Introduction 

Since the outstanding result of Montroll, Potts, Ward [1] there appeared many papers de- 
voted to the problem of the spin-spin correlation function in the two dimensional Ising model 
(IM) with the nearest-neighbour interaction on the infinite size lattice (see e.g. refs. in [2]). 
The achievements in this field are mainly connected with the analysis of the scaling limit 
[3], [4], [5], because IM just in this limit is of interest from the quantum field theory point 
of view. The so-called form factor representation [6] for the two dimensional IM correlation 
function appears to be a crown of this activity: 
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In the ferromagnetic region of the coupling parameter the summation in (1.1) is extended 
over even n, in the paramagnetic — over odd. It is worth of noting that this representation 
was first evaluated [6] in the framework of the S'-matrix approach [7] and then [5] by means 
of straightforward IM solution. The discovery of the form factor representation was very 
fruitful for the progress of exactly integrable quantum field theories [8]. 

The advantage of representation (1.1) consists in that the dynamical aspects of the system 
are separated from the kinematical ones: form factors squared are integrated over the phase 
volumes of the n-particle intermediate states. The form factor representation (1.1) clarifies 
the dynamics of the model but does not answer on the questions about the spectrum, sort 
and statistics of the particles that form the intermediate and asymptotic states of the system 
under consideration. The analysis of the model at finite (nonzero) temperature or in the 
finite volume is necessary for this purpose. One can observe an activity in this field in last 
years [9]-[15]. These works show in particular that the problem is complicated enough: the 
authors of refs. [14], [15], for example, call the conjectures of [11], [13] in question. I think 
that a simple exactly solvable lattice model example would be very useful for the business. 
Formally the finite temperature (in quantum field theory sense) means the finite size along 
the temporal axis and periodical boundary condition for boson fields or antiperiodical — for 
fermion fields. Meanwhile there is no representation which is analogous to the form factor 
one for the correlation function on the finite size lattice even for the Ising model. 

This work, I hope, makes up a deficiency. I have calculated the IM spin-spin correlator 
on the lattice wrapped on a cylinder by the use of the classic methods of IM theory [16] 
adapted properly to the case of the finite sizes. The solution is expressed in the form similar 
to the form factor expansion (1.1)— (1.3). For reader's convenience I write the result just in 
the Introduction. If one considers it obvious he would be free of cumbersome mathematical 



transformations. So, the expression obtained in this work is 

(a(r 1 )a(r 2 ))=tt T e- r / A Y,9n(r), (1.4) 

n 

*w=5wfny* (i.5) 

F u - TT sin( ^' ~ ft )^ 2 ) n fi^ 

More precisely, the representation (1.4)-(1.6) corresponds to the correlation of spins posed 
on the line parallel to the axis of a cylinder with N sites on its circumference. The values 
£, £ T , A, 7j, ?7j, in eqs. (1.4)-(1.6) are defined further in Sections 2, 3, 4. Three of them - 
£t, A and ?7j -- are specific cylinder circumference dependent values: ln£r, A -1 , ^ vanish 
if N tends to infinity. The appearance of the summation in (1.5) instead of the integration 
over phase volume is natural consequence of the size finiteness. It is widely believed that the 
underlying IM field is the fermion one: for instance the IM partition function is exactly the 
same as in the free Majorana fermion system. So, the boson spectrum of quasimomentum in 
Brillouin zone (this is denoted by the upper index in sum (1.5)) is surprising . Nevertheless 
it follows unambiguously from the calculations. It takes attention also that the lattice form 
factor (1.6) does not depends explicitly on the cylinder circumference. Note that similar 
expression for F n [q] at even n was found by authors of [5] for the infinite lattice case. 

In Section 2 the model is formulated and the brief evaluation of the Toeplitz determinant 
representation for the IM correlation function is given. The corresponding Toeplitz matrix 
allowing for the finite size lattice is calculated. In Section 3 the lattice form factor represen- 
tation for the correlation function is derived for the ferromagnetic domain of the coupling 
parameter and for the paramagnetic domain — in Section 4. The scaling limit is evaluated 
in Section 5. In Conclusion I discuss the relationship between the IM correlation function 
on a cylinder and quantum field Green function at finite temperature or in finite volume. I 
also comment on some prospect of generalization of the obtained results. In Appendix the 
factorised representations for the Toeplitz determinants of special form which are used for 
obtaining lattice form factor expansions are evaluated. 

2 The Model 

The Ising model with the nearest neighbour interaction on the square M x N lattice (see. 
Fig. 1) is defined by the hamiltonian H[a] 

-(3H[a] =lC£a(r)(V x + V>(r), K = (3 J, 

r 

where J is the coupling parameter, j3 is the inverse temperature; two dimensional vector 
r = (x,y) numerates the lattice sites: x — 1, 2, . . . , M, y = 1, 2, . . . , N; V x , V y are the 
one step shift operators 

V x <r(a;, y) = a(x + l,y), V y a(x, y) = cr(x, y + 1). 
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Fig. 1. The positions of the correlating Ising spins on the lattice. 



The partition function of the model is 



* = E' 



-f3H[a\ 



the pair correlation function is 



(2-1) 



(2.2) 



For the lattice with the periodical boundary conditions (wrapped on a torus) the par- 
tition function (2.1) can be expressed [17], [18] through the combination of the four Gauss 
Grassmann functional integrals with different boundary conditions 



Z = (2 cosh 2 K) MN UQ^ + Q(V) + Q W>) - g(»); 



where 



Q= d[^}e sw , d[t/j] = n# 1 (r)# 2 (r)dV 3 (r)# 4 ( r )- 
J r 

The action S[ip] is the anticommuting quadratic form 
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t = tanhK. (2.5) 
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The indices (/, b) in eq. (2.3) shows the type of the boundary conditions for the shift 
operators in the matrix (2.5). 

7" - antiperiodical: (V^) M = -1, (Vf») w = -1; 

the corresponding quasimomentum runs over halfinteger values in Brillouin zone 
(—71, 7r) (in 2n/M unites for p x and 2n/N - for p y ) 

"6" - periodical: (V^) M = 1, (Vff = 1; 

the corresponding components of quasimomentum runs over integer values. 

The correlation function (2.2) can be also presented in terms of the Grassmann functional 
integrals with the Gauss distribution [18] 

si/" + si" 1 + s!/" - ef' 

where Qd is the functional integral 

Q d = [dty}e s ^. (2.7) 



The action with the defect denoted by Sd[ip] in (2.7) differs from (2.4) in that the parameter 
t is replaced by t _1 on the links along the path which connects the sites i^ and r 2 . 

I consider the special case when the sites I - ! and r 2 are situated on the line parallel to 
the horizontal axis as is shown on Fig. 1: 

r 1 = (x,y), r 2 = (x + r,y), \r ± - r 2 | = r. 

The action with the defect in this case has the form 

sM = sm + Asm, 

where 

r 

AS[ip] = (r 1 -t)J2^\ x + x '- 1, y)^P 2 (x + x', y). 

x'=l 

The ratio in the r.h.s. of (2.6) simplifies 

Q (ff) 
(a^air,)) = f ^[1 + o(e~ M ' N )i (2.8) 

if at least one of the lattice sizes is much larger than other, i.e. torus degenerates to cylinder: 
in our case M/N ^> 1. Note, that only the terms with the antiperiodical boundary conditions 
survive. 



After some manipulations with the functional integrals the ratio in the r.h.s. of (2.8) can 
be transformed to the functional integral over Grassmann field which is defined just on the 
line connecting the sites ri and r 2 

((T(r 1 )<r(r 2 )> = f d[ x }e^ A ^ = det A, (2.9) 

r-1 

d ix] = II d xl dxl, ix^X 2 ) = J2 xlA x , x al>, 



x=0 



i v (/) eg-^!2Ml + t2 ) - ( x - t2 )( e ^ + t2 *~ ipx )] (9 im 

x ' x ' MN 2 ^ (1 + £ 2 ) 2 - 2t(l -£ 2 ) (cos p x + cos p y ) ' l j 

x, x' = 0, 1, . . . , r — 1. 

I recall that the index (/) points the antiperiodicity: the momentum's components run over 
halfinteger values - fermion spectrum 

E (/ V) = X>(|(' + !))> E (/) ^,) = X>(f (' + 1 

Px 1 = 1 Py 1 = 1 

Now let us rewrite the expression (2.10) in appropriate for the Toeplitz determinant 
theory manner. From the tabulated formulas [20] 

n (/) 2(cosh 7- cos Py ) =e^(l + e- JVr ) 2 , ^^(cosh^ -cosp y ) = e N ~<(l-e- N ~<) 2 , (2.11) 

Py Py 

it follows in particular 

1 y(J) 1 = tanh(JV 7 /2) 1 ^(6) 1 = coth(JV 7 /2) 12 

A^p cosh 7 — cos p y sinli7 iVr^ cosh 7 — cos p^ sinh7 

At the cylinder case M — > 00, iV = const the summation over p x turns into the integral 

i-£ (/) „(e--) = ± I *1 + (e-^). (2.13) 

M „ 27H J Z 

Px 1*1=1 

Summing up (2.10) over p y by use of (2.12) and accounting for (2.13) one can obtain for the 
matrix A XjX r 

A x , x , = J-<f -z-*+ x 'A(z), (2.14) 

where the kernel A(z) is 

A(z) 
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(1 — e~f*z)(l — e^^z^ 1 ) 
v 7 (I 1 -T(z), z = e~v, (2.15) 



and 

T(z) = tanh[(A^ 7 (p))/2], (2.16) 

70 = \2K + \nt\, -f n = 2K-\nt. (2.17) 

The function ^(p) in (2.16) is denned by the following equation 

4sinh 2 ^ = // 2 H-4sin 2 |, (2 .i8) 

where 

/j, 2 = 2(v / sinh2K - l/v / sinh2K) 2 . (2.19) 

It follows from eqs. (2.18), (2.19) in particular 

7(0) =70, i(^)=1tt, cosh7 7r -cosh7 = 2. 

The function T(z) (2.16) tends to unity at the thermodynamic limit N — > 00 (more 
precise Nfj, ^> 1), and the kernel A(z) (2.15) turns into the classic IM theory expression. In 
this limit the value K = K c is the critical point: the specific heat and correlation length 
diverge. At the finite iV the corresponding singularities are smoothed and there is no phase 
transition. Nevertheless also in this case the value K = K c is notable. The calculation 
technique based on Wiener-Hopf sum equation "distinguishes" the ferro- and paramagnetic 
regions and the expression for the correlation function at K > K c differs from that at 
K < K c . Just in this sense I shall mean further the words "critical point", "phase" and so 
on. Note that the following equalities take place at the critical point 

sinh2K c = l, /! = 0, 70 = 0, 7* = 21n("\/2 + 1). 

3 The Ferromagnetic region: K > K c 

The dependence of the correlation function (2.9) 

(<7(ri)<r(r 2 )) = detA< r ) (3.1) 

on the distance r = |ri — r2| enters through the matrix dimension denoted by upper index 
(r) in the r.h.s. of (3.1). The dependence on the cylinder circumference N comes into the 
function T(z) (2.16) and the dependence on the coupling parameter K is defined by the 
parameter \x (2.19) in the eq. (2.18) for the function 7(7?). 

The kernel (2.14) of the matrix (2.9) A^ has to be represented in the factorized form 
according to the method of the determinant calculation which we use (see (A6), (A7) in the 
Appendix), 

A(z) = P(z)Q(z- 1 ), 

where the functions P(z) and Q(z) are analytic ones inside the \z\ — 1 circle. It can be 
formally done by use of the projection operators (A8) that are defined in the Appendix 

2\n P(z) = \n(l - e~^ z) - \n(l - e^° z) +2V In T(z), 
21nQ(^ 1 ) = ln(l-e" 70 ^ 1 )-ln(l-e^ 7 ^" 1 ) + 2PlnT(2). 



Otherwise accounting for the explicit form (2.16) one can express T(z) through the ratio of 
the products with help of eqs. (2.11) 

II [cosh 7(p) — cos q] 

T 2 (z) = -i^ . (3.2) 

II [cosh 7(2?) — cosg] 
1 

It follows from (2.18) 



K 2 



cosh7(p) — cosg = 2 H cosq — cosp = ^^'(l — e 7 ^z)(l — e J ^ q 'z 1 ), 

so eq. (3.2) turns into the form 

n W e 7(9)(i _ e -T(9)«)(l - e~^ q h~ l ) 
T ^' = n (/) e^)(l - e-7(8)z)(l - e-^z" 1 ) ' 

Owing to the momentum g runs the entire Brillouin zone all the factors of the r.h.s. (3.3) 
are double except for the values q = and q = it in the nominator. So the functions A(z), 
P(z) and Q(z~ l ) have no branch points rather the simple poles and zeroes: 

n W( e 7(9) _ z ) n (6)(! _ e"7W z -l) 

D / \ /70 — 7vr\ 0<<?<7r „, _!s 0<g<7r , Q .v 

P W = Bq ,^_ r -j n , /1(e7 ,„_, ) . «* )= nM(1 _ e -,„ lrl ) < 3 - 4 > 

0<g<7T 0<g<?r 

According to (A25), (A28) 

CO 

detA«=e^X>,(r). (3.5) 

Let us write the function h(r) as the sum of three terms 

h(r) = -r/A + In £ + In £ T . (3.6) 

It follows from the eqs. (A23) and (2.15) 

1 rdz 

\ J z 



A" 1 = -[lnP(0) + hiQ(0)] = -^- <f — \nT(z) = - [ dp \ncoth(N 1 (p)/2). (3.7) 

2m J z n J 



The parameter A has the sense of the coherence length. Its asymptotic behaviour at large 

A^ is 

A^7o/__Z^. (3 . 8) 

V2sinh7o V ; 



One can see from (3.8) that the coherence length grows very rapidly when the cylinder 
circumference increases and 70 7^ 0. The lattice can be considered as infinite at Nyo ^> 
max[ln(iV7o),ln(r7o)]. It follows also from (A23) that 



H^ T ) = ^-Jdz\nQ(z~ 1 )^-\nP( 
or separately for the In £ and In £r 
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|Z2|>|*l| 



[Z2 - Z lj 



(3.9) 



ln£ = - In 
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In 6 



sinh 70 • sinh 7^ 

sinh 2 (( 7o + 7.)/2) 

d(z 1 z 2 



(2m) 2 I (z 2 - z x ) 2 

\z2\>\zi\ 



\nT( Zl )\nT(z 2 ). 



(3.10) 
(3.11) 



The r.h.s. of the eq. (3.10) does not depend on N and is singular at the critical point 

ln£~-ln/* at fi->0. (3.12) 

For the calculation ln£r let us return to the momentum variables and integrate (3.11) by 
parts. The result is 



or 



In 6 



In 6 



N 2 



dpdqi(p)i(q) 



2n 2 J sinh(iV7(p))sinh(A^7(g)) 



N 2 



d"f p d'y, 



pi* lq 



2tt 2 J smh(Ny(p)) sinh(N>y(q)) 

70 



In 



In 



sin((p + q)/2) 



sin((p-g)/2) 
sin((p + q)/2) 



sin((p-g)/2) 



Differentiating eq. (3.13) with respect to \x one can obtain 



d 1 

— ln£ T = --[cl(fi) - c 2 2 (fi)], 



where 



Ci(/i) 



N 



d^ / cosh 7^ — cosh 7 



2n J sinh(iV7) y cosh 7 — cosh 70 



(3.13) 



(3.14) 



(3.15) 



70 



N r 
CM = 2^S 



dy / cosh 7 — cosh 7 



70 



sinh(iV7) y cosh 7^ — cosh 7 



Deriving eqs. (3.14), (3.15) one has to account for that the integrand of (3.13) vanishes on 
the boundaries of the integration region. Note also that the derivatives of p and q are 



dp 



dq 



d/j, 2 2sinp' dfi 2 2 sing' 

9 



that follows from (2.18). For N'-fo ^> 1 the asymptotics of ci(/j,) and C2(/x) are 



Cl( ^"cosh(W)V2^h^' C2(/i) "4cosh(iV7o)V^- (5U6) 

Collecting (3.14)-(3.16) one can obtain 

00 
lnfr = \ dvv\c\{y) -c\{v)\ ~ -e" 2 ^ . 

J IX 

So, we see that £ T — > 1 if iV — ► 00. 

The integral (3.13) diverges at the critical point. This means that ln£r as a function of 
/j, has a singularity. It easy to obtain at // — > 

<91n£ T 1 

Cl(/i) " 1/2 ^' ^^ = -8y 

and consequently 

ln£ T = — - In (jb + const. (3-17) 

Therefore, one can see from (3.12), (3.17) that the function ln(££r) has no singularity at 
H — if N is finite. 

Let now evaluate the expansion in the r.h.s. of (3.5) 

00 
G F (r) = Y.92i(r). (3.18) 

1=0 

The coefficients g2i{r) are expressed through the 2/-multiple contour integrals (A30) as it is 
shown in the Appendix 

21 1-1 1 1 

r_iv r n(^on n [(^-i - ^_i) 2 (^ - Z2j ) 2 ] n^(^-i) 

£w) = ,,,,/o ^9/ f ; — i i • ( 3 - 19 ) 



m{2m)2l J<i n n (1 - ^-i^) 2 n w^ 1 ) 

j=l jr' = l 1=1 

The integral (3.19) is defined by the concrete form of the function W(z) = P(z)/Q(z^ 1 ): in 
our case this is the eq. (3.4). The integrand in the r.h.s. of (3.19) is analytic at \zi\ < 1 
except for the singularities which the functions W(zi) and W~ 1 (z~ 1 ) possess. To separate 
explicitly the singularities inside the circle \z\ — 1 these functions ought to be rewritten in 
the following way 



W{z) 



W'Hz 



i/~-i x 



P(z) 


"\ 


(1 - e~i"z)(l - e-^z- 1 ) P$(z) 


Q(z- 1 ) 


( 


1 - e-T°z)(l - e-^z- 1 ) T(z) ' 


Q(z) 




(1 - e-r°z)(l - e-^z- 1 ) Q 2 T (z) 



P(z- 1 ) \{l- e~^z)(l - e-^z- 1 ) T(z) ' 
10 



(3.20) 
(3.21) 



The analytic at \z\ < 1 functions Pr(z) and Qt(z) are defined by the equations 

f r dz' 

\nP T (z) =V\nT(z) = — A- lntanh(A^ 7 (p')/2), (3.22) 

2m J z — z 

kl=i 

1 r dz' 

\nQ T (z- r ) =V\nT{z) = — A -rntanh(7V 7 ( p ')/2), (3.23) 

2m J z — z 

\z\=l 

z' = e~ ip '. 
Note that it follows from (3.7), (3.22), (3.23) 

lnP T (0) = -1/A, lnQ T (0) = 0, P T {z) = P T (0) ■ Q T (z). 

With allowance of these eqs. one can obtain instead of (3.20), (3.21) 

W (z) = e^-^/ 2 - 2 / A (cosh 77r - cosp)Q 2 T (z) COtH ^ iP \ /2) , (3.24) 

sinh 7 (p) 

W-\z- 1 ) = e^-^/ 2 (cosh 7o - cosp)g|(^) C ° th( f 7(p) / 2) , (3.25) 

smh 7 (p) 

coth(JV 7 (p)/2) _ 1 y( b ) 1 

sinh 7 (p) Ng cosh 7 (g) — cosp 

It is seen from (3.24)-(3.26) that the singularities inside the integration contours of (3.19) 
are the simple poles at the points 

27r/r 

corresponding to the boson spectrum of the quasimomentum q. Accounting for 

1 rdz v(z) v(e _7(9) ) 



2m J z cosh 7 (g) — cos p sinh 7 (g)' 

|z|=l 

where v(z) is analytic at \z\ < 1, the integral (3.19) is fulfilled by residua. The result is 

vM-t^^ni-^-^M (3.27) 

where 

i=l j=i+l 

I 

I! [(1 + cosg 2i _i)(l - cosg 2i )] 
x ^ — » ( 3 - 28 ) 

n eT* n n (i - e^-™) 2 

i=l i=l j'=l 

11 



7* = 7(?i), ^ = ^?i) = -lBQr(e- 7(ft) ), (3.29) 



[g] 9i 92 <32i 



Cf = ,,,, — is the binomial coefficient. 



Using the equalities 



|(e 7l — e 7i )(l — e 7i+7j ) = cosli7j — cosh7j = cosg^ — cos<ji, 
the eq. (3.28) can be transformed to 

yzi[q\ = — t , 

nsinh 2 (( 7j + 7j )/2) 

Kj 

where 

i~i i 
V 2 i[q] = 2~ 21 Cf l Y[ II [(cosq 2i -i - cos q 2j ^ (cos q 2i - cos q 2j ) 2 ] x 

i=l j=i+l 
I 

xn(l + cosg 2l _ 1 )(l-cosg 2i ). (3.30) 

i=i 

The expression (3.30) is not symmetric with respect to changes the summation variables ^ 
with even subscripts by that with odd subscripts. The summation extracts the symmetric 
part of the function V 2 i[q] which is denoted by V 2 s / m [g] 

u i p 

where the sum over all permutations of q 2 i-\ and q 2 j (theirs number is just Cf 1 ) is denoted 
by J2 ■ The equality 

V 2 T[q] = U^riQ], (3.31) 

is the crucial for the final result. Here the even (with respect to qi <-> —qi) part of the 
product 

U a [q] = f[_8m 2 ( <k - qj 

i<j 

is denoted by U^^lq] 

[cr=±l] i<j 



12 



Consequently the function V2i[q] under summation (3.27) can be replaced by the function 
U2i[q]- This makes it possible to present the coefficient g2i(r) in the form similar to that of 
eq. (1.2) 

n = 2l, 

9 »M = ^rE W n(- E¥ )^fe]. (3^32) 

where F n [qi] is the lattice analog of the form factor (1.3) 

sm((q i -q j )/2) 



Fnh) = n 

i<j 



sinh(( 7i + 7i )/2) 



, F [q] = l. (3.33) 



I emphasize once more that the momentum spectrum over which the summation (3.32) is 
extended (the intermediate states summing up) occurs to be boson one contrary to the 
fermion spectrum by which the initial values were defined, in particular the matrix A x x > 
(2.10). 

The function rji = i](qi) that is contained in (3.32), decreases rapidly when the cylinder 
circumference grows. It follows from its definition (3.29) 

1 /• dzlnT(z) 1 j dp( l-e~~ /(q) cos p) , W , T . ... . 

lb = : f 7fa) ~ / h H ^lncoth7V 7p /2. (3.34) 

irt J e 7W — z it J cosh 7(g) — cosp 

\z\=l 



g-A^o 



One can find at A^7 ^> 1 

^ q > ~ e i(q) - i ^/2tiN sinh 7o ' 
and 77(g) — > if N — > 00. 

4 The Paramagnetic region: K < K c 

The correlation function is defined uniform through the determinant det A^' in the whole 
region of parameters (including K < K c ) . But the method of its calculation is to be modified 
because the Wiener-Hopf sum equation technique, which is used, demands the Toeplitz 
matrix kernel of the appropriate form. Both the kernel and the logarithm of the kernel have 
to satisfy the Loran expansion condition. Meanwhile, the matrix kernel (2.10) is reforming 
compare to the ferromagnetic case (2.15), when K goes across K c 

A(z) = z- l B(z), 

where 

Biz) 



(1 -e-i°z)(l -e~^z) _, N , N 

V^), (4.1) 



\ (1 -e-To^- 1 )(l-e-T^- 1 ) 
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so, rather the function B(z) than A(z) possesses appropriate analytical properties. The 
matrix A^ now has the following form 

A x , x , = ~ i —z-^^Biz). (4.2) 

\z\=\ 

Factorizing the kernel (4.1) 

B(z)=P(z)-Q(z- 1 ), 

where the functions P(z) and Q(z) are analytic inside the circle \z\ < e 70 one can obtain 
P W = e*p( — j u(I){elM _ z y <H' )= n ,/) (1 _ e - 7M ,- ir «■»> 

0<g<?r 0<q<7r 

Comparing with the factorized representation in the ferromagnetic case one can see that 
the corresponding products (3.4) and (4.3) differ from each other by one term: the factor 
(e 70 —z) is appeared in P(z) and the factor (1 — e _70 ,2 _1 ) is disappeared in Q{z^ v ). It follows 
from (4.3) 

P(0) = P r (0) = e" 1 ^, Q(0) = Q T (0) = 1, 

where the functions Pt{z) and Qt(z) are the same as in the ferromagnetic case (3.22), (3.23). 
The matrix (4.2) has the structure (A31) considered in the Appendix. Therefore its 
determinant can be represented according to the eqs. (A32)-(A35) by the following way 

detA {r) = e h{r+1) F(r), (4.4) 

oo 

nr) = E/a+i(0- 

1=0 

The function h(r + 1) is given by eq. (A23). Writing it similarly to (3.6) and allowing for 
the difference between the definitions (4.3) and (3.4) for the corresponding functions P(z) 
and Q(z) one can obtain 

h(r + 1) = -{r + 1)/A + ln(^ T ) + ln(l - e" 70 " 7 ^), (4.5) 

where the values A, £ and £t are the same as in the ferromagnetic phase (3.7), (3.10) and 
(3.11). With the allowance (4.5) the eq. (4.4) can be rewritten as follows 

det A& = sinh(( 7o + ln )/2)e h ^G P {r), (4.6) 

where 

CO 

G P {r) = 2e- 1 / A -( 70+7 ^ 2 P(r) = ]>> m (r), (4.7) 

1=0 
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and for the coefficients g2i+i(r) in the expansion (4.7) one obtains from (A35) 

f-lV r 2l+1 

to+1 (r)=2e- 1 /A-(7o + 7.)/ 2 _l^I_ (27r , )2 m j n(^4)x 

|z»|<l % ~ 

n n (z2i-z 2j ) 2 u n (^_i-^_i) 2 u(z 2l w(z 2l )) 

X — . (4.8) 

The functions W(;z) = P(z)/Q(z~ 1 ) and W^~ 1 (z _1 ) = Q(z) / 'P(z _1 ) entering the integrals 
(4.8) have the following form in this case 

W(z) = 2e~ 2 / A ~^ + ^ 2 Q 2 T (z)smh 2 1 (p)^-y {b) - -^ , (4.9) 

w TW ' yFJ N^ cosli7(g) -cosp v y 

^V 1 ) = >°^V 2 Q^)l£ (6) — — i . (4.10) 

2 A/ g cosh7(g) — cosp 

It follows from (4.9), (4.10) that the integrand in (4.8) has only the simple poles inside the 
integration contour at the points 

2irk 
3 (fc)=e -T(«w), g (fc) = — , A; =1,2,..., iV; 

similar to the previous case. At these points, obviously, 

cos p = cosh 'j(q), cosh 7(p) = cos g, sinh 2 7(p) = — sin 2 q. 

So, the integration contours in (4.8) may be squeezed and the coefficient g 2 i+i{r) is expressed 
through the sum of residua 

e -<* + i)/A (t) 2j+y e -*-^ v [q] 



where 



%<j 



i-\ i 

\2 



V^ +1 [g] = Cf +1 2- 2 W I] II (cos fe -cosg 2 ,)~x 

j=l j'=i+l 
i H-l i 

x ii n ( cos ?2i-i - cosg 2 i-i) 2 n sin2 ^- ( 4 - 12 ) 

j=l j'=i+l i=l 

The expression (4.12) differs from the corresponding ferromagnetic one (3.30). But it occurs 
that also in the paramagnetic case the symmetric part of the function V 2 / + i [q] - with respect 
to change of each variable with even subscript q 2 i by each variable with odd subscript q 2 j-i 
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(the number of the permutations is Cf +1 ) - coincides with the even part of the function 
U2i+i[q] 

i<j v z 

The equality analogous to the eq. (3.31) takes place 
where 

v 2 s i?M = 7 ^TrY,V2i + M 

°Z P 

97-4-1 

[<r=±l] *<J V 

Consequently the function Va+i[g] can be replaced by the function U2i+i[q] under the sum- 
mation in the r.h.s. of eq. (4.11). We obtain eventually 

n = 2l + l, 
*M = ^E (6> n(^-)^M, (4-13) 

™ - um^y *«->• 

The functions 7« = 7(g«) and rji = rj(qi) are defined by the eq. (3.29). 

The form of the coefficients g n (r) and of the function F n [q] is surprisingly the same 
both the ferro- and paramagnetic cases. In the ferromagnetic case the correlation function 
expansion is extended over even n and for the paramagnetic case - over odd n 

oo 

(a( ri y(r 2 )> = ^^ T )e- r/A Y,92i(r), K > K c] (4.15) 

1=0 

(a( ri )a(r 2 )) = sinhp^)(^ T )e-/ A f>m(r), K < K c . (4.16) 

v z J 1=0 

Notice that the even values of the circumference N were assumed in the factorization 
procedure for the functions P(z) and Q(z~ l ) (see eqs. (3.4) and (4.3) ). For the odd N the 
Brillouin zone does not contain the q = tt point in the boson case. On the contrary this value 
appears in the fermion spectrum. It is not difficult to account for this detail, it is irrelevant 
for the final result. 
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5 The scaling limit 

The IM correlation length /i _1 diverges when the coupling parameter K tends to the critical 
value K c . The scaling limit means that both the number of sites on the cylinder circum- 
ference and the distance between the correlating spins tend to infinity provided that the 
corresponding scaling variables are finite: 

N — > oo, r — > oo, (x — > 0, 

fiN = v = const, fir = p = const. (5-1) 

Contrary to the limit iV — ► oo, \x = const the specific "cylinder quantities" ln£<r, A -1 , 
77(g) do not vanish in the limit (5.1). They survive and nontrivially depend on the scaling 
circumference v. 

Before the evaluation of these functions note that all summands in (4.13) contain the 
factors exp(— r r ){q)) which restrict the summation area by small values of momenta q <^ 1 
if r — ► 00. Consequently 



'V 2 , • 2 9 



7(g) = 2arsinh y — + sin - ~ y/i 2 + q 2 <C 1, 
and we have in the limit (5.1), for instance, 

I 7r p-n(q) oo „-pw(g) 

— V , - . . = V -— - + o(exp(-p///)), (5.2) 

where e is some positive constant < e < 1 and w(q) = \fq 2 + 1. The sum in r.h.s. of (5.2) 
means 

E * (?) = E «(27ta s /i/). 

9 k=~ oo 

The lattice form factor (4.14) in the limit q <C 1 is reduced to the continuous one (1.3). It 
does not depend explicitly on circumference v. 

The calculation of eq. (3.13) gives the scaling limit for ln£r 

, c y2 1 duj x duj 2 qi + q2 ,. Q , 

2tt z J smh{vuL>i) sinh(z/a; 2 J gi — q 2 

According to (3.14) and (3.15) one can obtain other representation for this quantity 

00 00 

lnfr = v 2 \dv'v'c\v% c[v) = - f .™ (5.4) 

J n J qsvah(uuj) 

The coherence length (3.7) in the scaling limit can be written in the following way 

A" 1 = fiA(v), (5.5) 
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where 

OO CO 

A(u) = - f dq In coth(z/cu(g)/2) = - / f, q - . (5.6) 

n J n J sinn i/w 

o 1 v ' 

And for the last quantity (3.34) which is specific for the cylinder we obtain the following 
scaling limit expression 

riM = -^ / 7 2 lnooth(i^(p)/2). (5.7) 

TV J u 2 (q) + p 2 

The two dimensional Ising model can be considered as a lattice regularization of some 
quantum field model in (1+1) dimensional (euclidian) space-time. In our case one of the 
dimension is infinite and other is finite with periodical boundary condition along it. One 
can define the corresponding renormalized two-point equal time Green function at finite 
temperature and in infinite volume 

G(p,u) = z' 1 (a(0)a(r)}. 

Here z is the wave function renormalization constant. We shall argue below that in our case 
this constant is 

z = 2£. (5.8) 

The connection between scaling variables p, v and physical ones is tuned by the equations 

v = m/3, p = mR, 

where m is the renormalized field excitation mass, f3 is the inverse temperature and R is the 
spatial distance between correlating fields. 

Collecting the corresponding formulas, we obtain the following form factor representation 
for the renormalized two-point equal time Green function at finite temperature 

G(p,v) = fr(v>~ pAiu) T,9n(p,v), 9o = l, (5.9) 

n 

9n(p,v) = ^£ (6) n(e-^-V^)m (5.10) 

F n [q) = ft^^, Fi = l, (5.11) 

r)i = r)(qi, is), Ui = uj(qi) = \Jq^ + 1, p = mR, v = m(3; 
n = 0, 2, 4, . . . for the ferromagnetic case, 
n — 1, 3, 5, . . . for the paramagnetic case. 
The functions £t(v), A(z/) and r)(q, v) are defined by eqs. (5.3)-(5.7). 
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It is easy to see that all "cylinder corrections" (5.3)-(5.7) exponentially decrease when the 
circumference v grows. The summation in (5.10) over phase volume turns into the integral 

if v — »■ oo 

l£ ( "-- Jdq 
i/V 2tt J H 

q — oo 

and the expressions (5.9)-(5.11) turn into the classic form factor representation (1.1)— (1.3) 
for the case of infinite continuous plane (quantum field at zero temperature). 

The form factor (5.11) is Lorentz invariant quantity. It can be rewritten explicitly in 
terms of Mandelstam variables 



F n [q] = U(l-4/ Sl 

i<3 



,1/2 



where 

Sij = (uJi + Ldj) 2 - (qi + qj) 2 . 

Accounting for this property it is easy to calculate the Fourier transformation of the Green 
function 

G(p 2 ) = f d 2 Re ipR G(p, oo) 

and to obtain the Lehmann representation for the propagator 

G{P 2 ) = lEUP 2 ), (5-12) 

dQ 2 
P 2 + Q 2 

A »( Q2 » - ^ /(&£)«<£;« - «>>(!>• - 4 (5 - 14) 



oo 



Up 2 ) = / -4^ A n(Q 2 ), Q 2 = Ql-Ql (5.13) 



(ran) 2 



-oo Xi=1 l ' Xi=1 X Xi=1 



For the paramagnetic case the summation in (5.12) is over odd n. In this case the function 
G(p 2 ) possesses the transparent structure of singularities in the p 2 complex plane [21]. The 
first term in the expansion (5.12) gi(p 2 ) has the simple pole at p 2 = —m 2 

9i(P 2 ) = -T—2- ( 5 - 15 ) 

With the allowance of (5.15) the following normalization of the propagator 

lira [G(p 2 )(p 2 + m 2 )] = l 

leads to the eq. (5.8) for the wave function renormalization constant. The next terms in 
(5.12) have the branch points at p 2 = —(nm) 2 , corresponding to the thresholds of the n- 
particles intermediate states. The presence of the simple pole as the lowest singularity in 
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the propagator is the necessary and sufficient condition for the existence of the asymptotic 
states corresponding to observed particles. 

The next important question is about the sort of the particles: more precisely, about the 
type of statistics (Bose or Fermi) of the particles in the asymptotic and intermediate states. 
To clarify the problem we have to compare the IM result with that for the free boson and/or 
free fermion fields. Let us recall the expressions for the pressure and energy density of the 
free boson and fermion relativistic gas in (1 + 1) dimension 

p(/> fe ) 1 7 dug £^ ijdquo 



m? -k J e UU) ± 1 ' m 2 n J e™ ± 1 

1 o 

where v = j3fi, uj(q) = y/q 2 + 1; the signs (+, — ) correspond to fermion and boson respec- 
tively. The corresponding IM thermodynamic quantities are exactly the same as for the free 

fermion gas, i.e. 

p(i) = P U) £(I) = £(/)_ 

This was one of the reasons to interpret the IM as the free fermion system. 

By the way, one can see that the specific "cylinder quantities" (5.4) and (5.6) may be 
expressed through the free Bose- and Fermi-gas thermodynamic characteristics 

A(i/) = — 2 (P^ + P^), 

m 

The relativistic gas is the system with the zero chemical potential: the number of particles 
is not fixed rather it is defined by the thermodynamic equilibrium condition. As the conse- 
quence, the difference between fermion and boson gas at low temperature is quantitatively 
negligible because the gas is dilute and becomes apparent only at high temperature. For 
example, 

P {b) - P (f) f e-72 3 / 2 for v > 1 

P (b) + P (f) ~ J 1/3 f or v <^\. t 5 - ib J 

The free field propagator & b '^\p 2 ) is the simple pole (is proportional to the simple pole for 
the fermions) 

&^(p 2 ) = ^- 2 , P 2 =pl+P 2 v 

At finite temperature the zero component of the momentum becomes discrete 

(fc) _ J 2irk/(3 for bosons 

Po = \ 2rr(k + l/2)//3 for fermions 

In result the spatial correlations of the free boson and fermion fields are 

&J> M =•£*»££. (5.17) 

q 2vuj{q) 
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The corresponding IM pole contribution is 

G%( P , v) = frfr)e-^X; ( > ) e 2 ^ (g) ■ (5.18) 

At low temperature z/ ^> 1 the functions ln£r(z/), A(z/) and r](q, v) are exponentially small 



v 



c(v) ~ e- u J2/7ris, ln£ T (z/) ~ — e"^, 

V 7T 



A(z/) ~ e^y^W, V {q,u) 



4e" 



cj(g)v2vn^ 



and the function Gp 0le (p, f) coincides with that for the free boson field. 

The high temperature behaviour is more instructive, as was noted above. This assertion, 
at first sight, contradicts to our nonrelativistic quantum physics experience: the lower tem- 
perature - the higher role of the quantum effects. In fact, the quantum corrections become 
apparent if the occupancy of the energy levels becomes large enough. If the particle density 
is fixed this is realized when the temperature tends to zero. But in the relativistic gas the 
particle density decreases exponentially, due to the annihilation, when the temperature di- 
minishes. Hence, the occupation numbers tend to zero and the difference between the Bose 
and Fermi statistics is irrelevant. On the contrary, at high temperature the level occupancy 
depends on the competition between the extension of the accessible levels and the growing 
particle density. The quantitative issue of the described picture is given by the eq.(5.16) for 
the pressures. The same is true for the correlation functions. One can see from eq. (5.17) 
that for v >> 1 (low temperature) the difference between boson and fermion correlation func- 
tions is quantitatively negligible. But for v <C 1 (high temperature) the difference between 
boson and fermion field correlations is appreciable, as it is seen from the following equations 

G (b) (p,u) ~ — e~ p + — e~ 27r ^, 
2u 2n 

G (f \p,u) ~ -*-*pI v . 

71 

The high temperature asymptotic behaviour of the functions £,t{ v ) is 

/ 2 \V4 

\7lW 

where 

oo 

C = expf — f —e~ 2x t&nh 2 x\ ~ 1.04. (5.19) 

o 
Note that the integral in (5.19) can be expressed through the Glaisher's constant Cq so that 

C = V^e 1/4 2^ 6 C G 3 . 
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The asymptotics of A(u) and r)(q, v) for v <C 1 are 

A(v) ~ — e-**") - ^^ 

In result we obtain from (5.18) 

G%( P ,u)^C^) 1/4 e-^^. (5.20) 

One can see that the IM correlations (5.20) decrease slower than the free fermion gas corre- 
lations, if the temperature grows, but more fast (the screening) compare to the free boson 
gas correlations. 

Therefore, on the level of the finite temperature Green function the quantum field system, 
corresponding to the IM scaling limit from the paramagnetic region, looks like strongly 
interacted bosons. The thermal fluctuations dress the intermediate states: the everaging over 
phase volume is accompanied by the additional thermal weight function exp(— r](q, v)). But 
the dynamic quantities (form factors squared) are left bare. In spite of this the corresponding 
thermodynamic quantities (free energy, pressure, energy density) are identical to those of 
the free fermion gas at any temperature - the miracle in a sense. 

The picture of the ferromagnetic region scaling limit is more sophisticated. First of all the 
nonvanishing vacuum expectation value at zero temperature appears (a) ^ 0. So, the Green 
function of the excitations over the condensate is to be defined. In our case it means that 
the first term in the expansions (5.9), (5.12) is to be removed and summation is expanded 
over n = 2,4, . . .. The lowest singularity of the propagator is contained in 52 (p 2 ) term of 
the Lehmann representation (5.12)-(5.14). It is not simple pole. Consequently there are no 
asymptotic states corresponding to the particles. From the other hand the thermodynamic 
quantities here are the same as in the paramagnetic case: the matter exists, the particles do 
not. This paradox is spurious. Due to the self dual properties of the IM there exists side by 
side the order parameter (a) the disorder parameter (/i) . The Green function corresponding 
to the disorder parameter in the ferromagnetic region is identical to that corresponding to 
the order parameter in the paramagnetic region. So, there exist the quantum field excitations 
corresponding to the particles states also in the ferromagnetic region. But these states cannot 
be generated by the external sources which are locally connected with the quantized field 
a. Moreover, the corresponding sources apparently do not commute with the internal field. 
The picture seems to be instructive, in particular, for QCD confinement problem. 

Note also that the correlation function on a cylinder in the ferromagnetic region K > K c 
vanishes if p — ► 00 at any finite circumference v 

lim G(p, v) = { J £ n 

p^oo Kl ' ' y for v > 

Therefore, the IM field condensate appearing in the ferromagnetic region is not stable. Ther- 
mal fluctuations at any low temperature destroy them. 
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6 Conclusion 

It should be observed first of all that, although the representations (4.15), (4.16) for the 
correlation function are obtained for different regions of the coupling parameter (K > K c 
and K < K c respectively) , both of them are valid and equal to each other for all values of 
K at any finite cylinder circumference N. For the finite N the expansions (4.15), (4.16) are 
not series rather finite sums due to the form factor F n [q] is identically zero for n > N. Only 
the thermodynamic limit N — »■ oo splits the domain of their validity. 

The obtained lattice form factor representations possess one beautiful feature: they can 
be checked numerically or analytically. For small values of N the result may be computed 
in other way, for example with the help of the transfer matrix technique. In particular at 
N = 1 we have the one dimensional Ising chain. The elementary calculation gives 

(a(0)a(r))=f. (6.1) 

One can see after simple transformations that both the representations (4.15) and (4.16) 
coincide with (6.1). 

Apart from its intrinsic interest the Ising model deserves scrutiny because of its relevance 
to nonperturbative explorations in quantum field theory. Mutual correspondence between 
field model and spin lattice system in scaling regime implies in our case the limit K — > K c , 
(N,M,r) -^ oo at {Nfi, Mfi, rfi) = const. Mysterious process of thermodynamic or 
scaling limit may be observed in detail proceeding from the explicit expressions obtained in 
this work. 

The IM correlation function (a(r)a(r')) in scaling limit corresponds to the quantum 
field model two point Green function G(x,t; L,/3). The variables of the Green function are 
assumed to be scaled on mass: m = 1, x = fi\r x — r' x \ is the spatial distance, t — n\r y — r' \ 
is the temporal one. L = \xM is the volume, (3 = fiN is the temperature. The boundary 
conditions along the temporal coordinate are quite definite: periodic for the boson field, 
antiperiodic for the fermion field. This is the issue of the appropriate formulation of the 
quantum field theory at finite temperature. On the contrary, the spatial boundary condition 
may be arbitrary. In any case the presence of boundaries breaks down explicitly the Lorentz 
invariance (the rotational invariance in Euclidian metric case). The spatial and temporal 
coordinates are not equivalent in general. In this sense the configuration of correlating 
spins considered in this paper corresponds to the equal time Green function G(x, 0; oo,/3) 
in the infinite volume at finite temperature. With some stipulations this function may be 
viewed also as that at zero temperature in finite volume L = /j,N (with periodicity in spatial 
coordinate). One has to keep in mind in addition that this is not equal time rather equal 
site Green function G(0,t;L, oo). It is necessary to consider the correlating spins placed on 
a cylinder circumference to compute the equal time Green function in finite volume and zero 
temperature. The corresponding problem may be solved by the method used in this work 
but with the Toeplitz determinant technique being properly modified. 

Meanwhile, the representation (4.13)-(4.16) for the IM correlation function has such 
transparent structure that the conjecture for the case of general displacement of correlating 
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spins suggests itself. The simplest one is 
(a(0)a(r x ,r y )) - 



Qnyxt Ty) 



(Z-tT)e- r * /A J29n(r x ,r y ), 



g-n/A^,^ _r>^ f e -r -,--,<■,,,-/,. 



n\N n 



y." n 



sinh 7j 



^M- 



(6.2) 
(6.3) 



It is amazing that this generalization is in exact agreement with the transfer matrix results 
for the N-rows Ising chains. We have checked it analytically for N = 2, 3, 4 and numerically 
for N = 5, 6. 

I thank V.Shadura for the fruitful discussions and O.Lisovy for the assistance in compu- 
tations. 
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Appendix 

The elements of the Toeplitz matrix arranged along the diagonals are equal 

A x , x > = M for k = x - x' = 0, ±1, . . . , ±r - 1 , (Al) 

\ 



A (r) 



( A 

A- L 
A 2 



A_i A_ 2 
A A_i 
Ai A 



A_ r+ i 

A_ r +2 
A-r+3 



(A2) 



y A r _i A r _ 2 A r _ 3 • • • A / 

The superscript in A^ shows the matrix dimension. If there is the function of the complex 
variable A(z) analytical in the ring a < \z\ < a -1 (a < 1) such that the coefficients A& of 
its Loran series 

oo 

A(z) = ]T A k z k (A3) 

k=— oo 

coincide with the matrix element (A2) then the matrix (Al) can be represented in the form 

A x , x , = ^-.i -z- x+x 'A{z). (A4) 

2m J z 

kl=i 

The function A(z) in the integrand of (A4) is called the kernel of the Toeplitz matrix A^'. 
If the inverse matrix is known then the calculation of the determinant reduces to calcu- 
lation of trace with help of the identity 



mdetA^(A) = /rfASp (A^(A)) _1 — A^(A) 



-id 



+ const, 
where A is a parameter of which the matrix elements depend on explicitly. 



(A5) 



24 



The inversion procedure for the Toeplitz matrix simplifies if both the kernel (A3) and 
In A(z) may be expanded in Loran series. If so, the kernel A(z) can be expressed in factorized 
form 

A(z) = P(z)-Q(z- 1 ). (A6) 

The functions P(z) and Q(z) are analytic inside the circle \z\ < a -1 : 



InP(z) = P\nA(z), lnQ(z _1 ) = P\nA(z). 



(A7) 



The projection operators P and P extract the analytic part of the given function v(z) inside 
and outside the circle \z\ — 1 correspondingly 



V 2 =V, V 2 = V, PP = PP = 0, P + P = l, 



Pv(z) 



2ni J t — z 2m J z — t 

\z\=l \z\=l 



(At 



The problem of the Toeplitz matrix inversion can be reduced to solving the system of 
the Wienner-Hopf sum equations [16]. In terms of the projection operators (A8) the matrix 
elements of the inverse matrix (A^)~ l x , are 



(^ (r) )^ = Vq^Pw-\z)[1 - V^z)Vw-\z)\- x V-^^ 



2=0 



(A^)-y = V^—Vw(z)[l-Vw-\z)Vw(z)r l V Z —- 



2=0 



(A< 



r)\-l 



P 

P(z 



4^y{l " V[l - w(z)Vw-\z)V}-Mz)r™-\z)}PQ^T 



2=0 



where 



w(z)=z r W(z), W{z) 



P(z) 



(A9) 
(A10) 
, (All) 

(A12) 



Each of the alternative but equivalent representations (A9)-(A11) for the inverse matrix 
{A^)~ l x , are preferable for the calculation of different matrix elements. In particular, the 
eq. (A9) is suitable for the calculation of the element (A^)~i 10 , the eq. (All) - for the 
element (A^)^: 



-lf-vM- 1 



(A {r) U = p (0 ! Q(0) {l " m - w(z)Pw-\z)P]- 1 w(z)Vw- 1 (z)} 



2 = 



(A13) 
• (A14) 



Expanding the inverse operators that are contained in the r.h.s. of (A13) and (A14) into 
geometrical progression 

oo 

[l-Pw{z)Pw- 1 (z)]~ 1 = Y\Pw{z)Pw- 1 {z)}\ 

1=0 
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[l-w(z)Pw- 1 (z)V]~ 1 = Y\w{z)Vw-\z)V] 1 



1=0 



and using the integral representations for the projection operators V, V (A8) one can obtain 
for (A13) with allowance (A12) 

1 oo 

(A15) 



(^ (r) );- 1 i,o = 


1 oo 

"Q 2 (o)S a2m(r) ' 




1 / 


21+1 

11 {dzizl) 

i=\ 


n w(z 2l ) 



^H-i(r) = ^-^r / ^ -P -■ (A16) 

n<i 2122+1 n (i - %+o n w(^i+i) 

i=l i=0 

In particular, the coefficients Oi(r) and 03 (r) are 

. 1 /■ dzz r ~ 2 

Oi(r) 



2m J W(V 



-IV 



Lz =1 



1 /■ d(z 1 z 2 z 3 ) (z 1 z 2 z 3 ) r W{z 2 ) 

o 3 (n - 



N<1 

Analogously one can obtain for the eq. (A14) 



(27Ti) 3 ( 7 Z!Z3(1 - Z lZ2 )(l ~ Z2Z3) W(Zi X )W(% X ) 



-1\ ' 



1 , 00 -. 

<- 4,r '»»"-?(0)O(0)( 1 -S a2 ' (r) )' (A17) 

2« Z 

° 2 ' (r) = TfcriV* / & *7 ( Al8 ) 

v ; i«ii<i Z1Z21 n (1 - 2i^+o n ^(<o 
»=i »=i 



and, in particular, 

1 /■ d(zi22)(«i^) r W^l) 



a 2 (r) 



a A (r) 



(2m) 2 7 ^ 2 (l-^ 2 ) W^ 1 )' 

F»|<1 

1 /■ g^i^^^)^!^^^)'' ^(21)1^(23) 



(2m)\ 7 2i2 4 (l - ^1^2) (1 - * 2 2 3 )(1 - ^4) ^(z 2 - 1 )^(^ 1 )' 

F»|<1 



Consider now the auxiliary Toeplitz matrix A^ r \\) depending on the parameter A, 

A x ^{\) = ^-. I - z - x+x 'P(z)Q x (z- 1 ), 0<A<1, (A19) 

lixi J z 

1*1=1 

where the functions P(z) and (5(2) are the same as in the kernel of the matrix A^. It is 
seen from the definition (A19) that at A = 1 
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A X ,A1) = ^~ f ^z- x+x 'P(z)Q(z- 1 ) = A x , x ,, 

and at A = 

A x , x ,(0) = 2^ / ^-z-* + *'P(z). (A20) 

|z|=l 

The matrix (A20) is the triangular one: all its elements to the right of the main diagonal 
are equal to zero 

A x , x (0) = P(0), A x , x ,(0) = if x' > x. 

Therefore 

detA^(0) = P r (0). (A21) 

The integration constant in (A5) is defined with help of (A21). 

Let us select the contribution h{r) in In det A^ which does not vanish when the matrix 
dimension grows 

lndetA (r) = h(r) + Ah(r), Ah(r) -+ if r -»■ oo. 

It is not difficult to observe that only the first term of the inverse operator expansion in (A10) 
contributes in h(r). This fact results in drastic simplification of corresponding computations. 
The inverse matrix (A10) is 



(-""•wo = vi^r^-v 2 



x r 



P{z) Q\z~ l ) P{z) 
1 r d(z 1 z 2 z 3 )z^ x {z 2 /z 3 )z 3 x ' P(z 2 ) 



(2tu) 3 J z 1 {z 2 -z l ){z 2 -z 3 ) P(z 1 )Q x (z 2 1 )P(z 3 )' 

\ z 2\>\zi,3\ 



9 A ^) = ^- f ^z-^'P(z)Q\z-i)\nQ(z-% 



\z =1 



The derivative of the matrix (A19) is 
d_ 

and finally one has for the trace in (A5) 
Sp 



(>«'"(A))-'^<"(A)] = pi^ / n>x 

P2i|>F2j'-l| 

x [1 - (z 4 / Zl ) r ][l - (z 3 /z 4 ) r }(z 2 /z 3 ) r P(z 2 )P(z 4 )Q x (z^)lnQ(z^) _ 
(z 2 - z 1 )(z 2 - z 3 )(z 4 - z 1 )(z 4 - z 3 ) P(zi)Q x (z 2 1 )P{z 3 ) 

= rlnQ(0) + — / dz\nQ( Z - l )^-\nP(z) + oir- 1 ). (A22) 

2m J oz 



z\=l 
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It is seen from (A22) that nonvanishing contribution in the trace occurs to be independent 
on the parameter A. Therefore, with allowance (A21) we obtain 

h(r) = r\\nP(0) + lnQ(0)] + — I dz\nQ( z' 1 )— In P(z). (A23) 

2iri J oz 

|z|=l 

For the calculation of Ah(r) the eq. (A5) may be used by keeping the next terms in the 
expansion of (A10). But it is more convenient for this purpose to exploit the specific property 
of the Toeplitz matrix i.e. 

M(r+D r i _ detAM 

1 Jo '° detAH- 1 )' 

Therefore 

detA (r) = detA (r+fc) f[ (A^)~J, 

s=r+l 

and after substitution (A17) we obtain 

r+k ,00 » 

detA^ = e h ^ r+k)+Ah{r+k) [P(0)Q(0)]- k J] (l-X>2i(s))- ( A24 ) 

s=r+l ^ Z=l ' 

It follows from (A23) that 

e h(r+k) = [p( )g(o)]V (r) , 

so, the factors -P(O) and (5(0) in (A24) cancel. Due to this fact the limit k — > 00 in (A24) 
can be taken 

detA (r) = G(r)-e h{r \ (A25) 

where 

OO • CO s 

G(r)= II (l-E«a(a))- (A26) 

Expanding the product of sums in the r.h.s. of (A26) 



OO/OO \ CO CO/CO 



G(r) = l- Y: E«a(*))+ E E E«a(*0 )(E«2i(*0 )-•••, (A27) 

si=r+l M=l / si=r+l s 2 =si+l V Z=1 ' V /=l / 

summing up over Sj and collecting the terms with the definite multiplicity of integration, 
one can represent (A27) in the form 

■CO 

G(r) = E0a(O- (A28) 

z=o 

The coefficients <?2z( r ) can be obtained from (A27) directly but it is more simple to compute 
them by use the following recursion relation 

i— 1 00 
0a(r) + E E <?2fc(s)a 2(i _ fe) (s) = 0. (A29) 

fc=0s=r+l 
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It is not difficult to deduce this relation from (A27). Using the eq. (A18) for a 2 i(s) and 
initial condition go(s) = 1 one can obtain from (A29) 



-V. 



92i{r) 



l\l\(2m) 



21 



21 l-\ I l 

n(^<)n n [(^-i - z 2j ^f{z 2i - z 2j y] uw(z 2i _ 1 ) 
J= i j=i+i i= i 



i=i 



in particular, 

92(r) = 

9±{r) = 



\zi\<l 



1 



n n(i-^-i^) 2 
i=ij=i 



n w(z£) 

4=1 



(A30) 



d{z t z 2 ) { Zl z 2 y W( Zl ) 

(27ri) 2 i a-z lZ2 y w(z 2 y 

\Zi\<l 

d(z l z 2 z 3 z 4 ) (z 1 z 2 z 3 z 4 ) r (z 1 - z 3 ) 2 (z 2 - z 4 ) 2 W{zi)W{z 3 ) 
A(2m)\ J (1- z x z 2 )\\- z 2 z 3 f(l- z 3 z 4 )i(l- z i z i yW{z 2 l )W{z 4 l \ 



Note that the integrand of (A30) is symmetric both with respect to permutations of the even 
subscript variables and with respect to permutations odd subscript variables. 

The structure of Toeplitz matrix is altered in paramagnetic region. Here it is necessary 
to evaluate the determinant of the matrix which has the following form 



B (r) 



( A_! A A x 
A- 2 A_! A 
A- 3 A_ 2 A_i 



B„ 



y A_ r A_ r+ i A_ r+2 
1 r dz 



4-2 \ 

4-3 
4-4 



(A31) 



2ni J z 
Ul=i 



x+x '- l A(z). 



The technique stated above is not valid because the kernel of the matrix B^ 

B(z) = z- l A{z) 

does not satisfy the Loran expansion condition for In B(z) if In A(z) does satisfy. Never- 
theless, the problem can be reduced to that considered above. Really, it follows from the 
definition of the inverse matrix element through the determinant and cofactor 

det£ (r) = (-l) r (A r+1 ))- 1 detA( r+1 ). 

Using the eqs. (A15), (A16) for the inverse matrix and the eqs. (A25), (A28), (A30) for the 
determinant one can obtain 



detfl (r) = (-l) r e h(r+1) F(r), 



(A32) 
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where 

1 oo 

F(r) = G(r + 1) ■ (A^ 1 ^ = ^yX>+i(r), (A33) 

i 
/a+i(r) = E a 2fe+i(r + 1)02(1-*) (r + 1). (A34) 

All summands in the r.h.s. of (A34) are the integrals of 11 + 1 multiplicity. After reduction 
the integrands to common denominator and symmetrization over even subscript variables 
and over odd subscript variables one can obtain 

/a+i(r) 



(A35) 



f-lV r 2l+1 

K J eh Y\ (r1~ ~ r ) -< 




/!(/ + 1)1(2™)^ f 11^-^- 

|z»|<l 

l-l I I l+l 

n n (z2i-z 2j ) 2 n n (^_i- 

j = l J=j + 1 4 = 1^ = 1 + 1 


-z 2j -i) 2 U(z 2l W(z 2t )) 

i=l 


X J / 

n n(i-^^+i) 2 

i = l 7=0 


U(z 2t+1 W(z^ +1 )) 

i=0 



and, in particular, 



, , , l r dzz r 
A( r ) = 7T- 



2ni J zW(z- 1 )' 



z =1 



, / x __ 1 / dfciz^s) (z 1 z 2 z 3 ) r (z 1 - z 3 ) 2 z 2 W(z 2 ) 



2(27Ti) 3 , / (l-«l22) 2 (l-^3) 2 (l-^l) 2 ^ 3 W(«r 1 )W(^' 1 )' 
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